NiO: Correlated Bandstructure of a Charge- Transfer Insulator 
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The bandstructure of the prototypical charge-transfer insulator NiO is computed by using a 
combination of an ab initio bandstructure method and the dynamical mean-field theory with a 
quantum Monte-Carlo impurity solver. Employing a Hamiltonian which includes both Ni-d and 
0-p orbitals we find excellent agreement with the energy bands determined from angle-resolved 
photoemission spectroscopy. This solves a long-standing problem in solid state theory. Most notably 
we obtain the low-energy Zhang-Rice bands with strongly k-dependent orbital character discussed 
previously in the context of low-energy model theories. 
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The quantitative explanation of the electronic struc- 
ture of transition metal oxides (TMOs) and other mate- 
rials with correlated electrons has been a long-standing 
challenge in condensed matter physics. While the basic 
concept explaining why materials such as NiO are insu- 
lators was formulated by Mott already a long time ago 
the development of an appropriate, material-specific 
computational scheme proved to be a formidable task. 
The electronic structure of the late TMOs, including the 
cuprate superconductors, is not only affected by the elec- 
tronic correlations, it is further complicated by the hy- 
bridization between the transition metal d-statcs and O 
p-bands located between the lower and upper Hubbard 
bands formed by the transition metal d orbitals. For such 
materials Zaanen, Sawatzky and Allen introduced the 
term " charge transfer insulator" , a prototypical example 
of which is NiO. In principle the simple crystal structure 
of NiO allows for a straightforward comparison between 
theory and experiment. However, a theoretical descrip- 
tion of the NiO bandstructure is made difficult by the 
competition between the local many-body effects, due 
to strong Coulomb interaction between Ni d electrons, 
and the band dispersion, due to the lattice periodic- 
ity, both observed with the angle-resolved photoemission 
spectroscopy (ARPES) 0,13. 

In this Letter we use a combination of a conventional 
bandstructure approach, based on the local density ap- 
proximation (LDA) , and the dynamical mean- field theory 
(DMFT) @,[6i,0| to investigate the bandstructure of NiO. 
No adjustable parameters enter. While the application 
of the LDA-FDMFT [1,11, [l^l framework has proven suc- 
cessful for the early TMOs, the charge-transfer materials 
were routinely avoided due to the additional complex- 
ity arising from the presence of p-bands. In the present 
work the O p-orbitals and their hybridization with Ni 
d-orbitals are explicitly included, thus allowing for a uni- 



fied description of the full spectrum. Our results reveal 
a non-trivial effect of the p — d hybridization in strongly 
correlated system studied so far only in terms of simple 
models [III, 111, [H . 

The application of the standard bandstructure theory 
to NiO is marked by a failure of LDA to produce an 
insulating groundstate [l3 |. The antiferromagnetic or- 
der within LDA [T^, despite rendering NiO an insulator, 
does not present much of an improvement since (i) the 
band gap is severely underestimated, (ii) the experimen- 
tally observed lower Hubbard band is completely missing, 
and (iii) static quantities such as the local magnetic mo- 
ment do not agree well with experiment. Moreover, the 
ARPES data of Tjernberg et al. [11] measured across the 
Neel temperature Tn—525 K show that the NiO band- 
structure is rather insensitive to the magnetic order. The 
first attempt to include the strong on-site correlations 
into the first-principles bandstructure methods was the 
LDA+U theory of Anisimov et al. 17|. The static, or- 
bitally dependent self-energy of LDA-I-U enforces a sepa- 
ration of the occupied and unoccupied d-bands and thus 
opens a gap of the experimentally observed size. This in 
turn leads to a significant improvement in the descrip- 
tion of groundstate properties such as the local moment 
or the lattice constant [l8| . Despite this success LDA+U 
does not provide a good description of the photoemission 
spectra since it places almost all of the valence d spectral 
weight into the lower Hubbard band. Exact diagonal- 
ization studies of Fujimori et al. (1^] on small clusters 
provide strong evidence that the dynamical correlations 
are necessary to capture properly the distribution of the 
d spectral weight within the valence band. 

A systematic inclusion of local dynamical correlations 
into lattice models was made possible by the dynami- 
cal mean-field theory (DMFT) @, H, |3|- The connec- 
tion of DMFT with bandstructure methods, usually re- 
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FIG. 1: (color online) The k-resolved total spectral function 
yl(k,aj) along the F-X (upper panel) and F-K (lower panel) 
lines in the Brillouin zone depicted as a contour plot. The 
symbols represent the experimental bands of Shen et al. [ij. 
The theoretical gap edge was aligned with the experimental 
one. 



ferred to as LDA+DMFT scheme 0,0, provides ac- 
cess to material-specific single-particle spectra as well as 
more general correlation functions. LDA-f DMFT calcu- 
lations usually start with the construction of an effective 
Hamiltonian from the converged LDA bandstructure us- 
ing a projection of the low-energy bands onto Wannier 
orbitals. Typically only the transition metal d bands 
are included in the construction, which amounts to in- 
tegrating out O p-states [19] at the price of having more 
extended Wannier orbitals, a procedure only justified for 
the early TMOs with low- lying p-bands. In the present 
study we project Wannier states out of a single energy 
window including both p and d bands, resulting in a more 
localized 0-p and Ni-d Wannier orbitals with mutual hy- 
bridization. 

The computation proceeds in several steps; (i) the 
construction of the effective Hamiltonian from converged 
LDA results, (ii) the self-consistent solution of the DMFT 
equations on the Matsubara contour, (iii) the analytic 
continuation of single-particle Green function to the real 
frequencies, (iv) the calculation of the self-energy on the 



real axis, and the computation of the generalized band- 
structure. We use a projection onto Wannier functions 
[20| to obtain the Hamiltonian discussed in Ref. 21. Next 
we solve the DMFT equations iteratively on the Matsub- 
ara contour, a key part of which is the auxiliary impu- 
rity problem handled by quantum Monte-Carlo (QMC) 
method [23 |. To ensure ergodic sampling we introduced 
global moves between ferro- and anti-ferromagnetic con- 
figurations of local Cg spins. The computational parame- 
ters are the same as in Ref. [2l|. The d-spectral functions 
were obtained by the maximum entropy method [23| ap- 
plied separately for the Cg and t2g symmetry. The local 
self-energy 'S{uj~^), which is formally an 8 x 8 matrix with 
the only non-zero elements on the diagonal of dd block, 
is obtained by solving the equation 



simultaneously for Cg and t2g symmetry. Here the Gdd is 
the diagonal element of the Green function corresponding 
to either Cg or t2g symmetry, is the 8x8 Hamiltonian 
matrix on a mesh of k-points and fi is the self-consistently 
determined chemical potential. The equation is solved 
approximately on the contour SJw"'' = O.ObeV subject 
to the constraints 3S(ci;+) < and Kramers-Kronig re- 
lations, which are, however, satisfied automatically for 
good quality QMC data. 

In Figs. [U we compare the theoretical bands, rep- 
resented by the k-dependent spectral density A(k,a;), 
along the T — X and T — K lines in the Brillouin zone 
with ARPES data of Ref. 4. Both theory and experiment 
exhibit two relatively flat bands at -2 and -4 eV followed 
by several dispersive bands in the -4 to -8 eV range and a 
broad incoherent peak around -10 eV. Overall we find an 
excellent agreement. The deviations around the F point 
in the lower panel of Fig. [T] are due to the inaccuracy in 
the location of the F point in off-normal-emission experi- 
ment . The crosses near the F point mark a weak band 
which was interpreted as a consequence of AFM order Q 
and is therefore not expected to be found in the param- 
agnetic phase investigated here. In Fig. [2] we show the 
orbital decomposition Ai„^(k,u!) of the spectral density 
visualized as a shape-preserving function 



order to capture both sharp and broader features in a 
single plot. 

We start the discussion of our results by considering 
the limit of vanishing p — d hybridization, in which case 
the entire d spectral weight is located in more or less fea- 
tureless Hubbard bands located below the p-band mani- 
fold [2lj and the holes in the uncorrelated p-bands have 
infinite lifetimes. The p — d hybridization changes this 
picture qualitatively. In particular, two additional bands 
of mixed character appear at -2 and -4 eV of Fig. [T] 
These bands contain about half of the valence d spec- 
tral weight as can be seen from the k-integrated spec- 
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FIG. 2: (color online) The orbitally decomposed spectral 
function A^i, (k, uj) along the F-X (left column) and F-K (right 



C+A„„{k,t^) ■ 



column) lines in the Brillouin zone plotted as 
The panels from top to bottom show the 0-p, Ni-d — eg and 
Ni-d — t2g contributions. Here C = 1.5, 2 for the p and d pro- 
jections, respectively. Detail of the uppermost valence band 
marked by the dotted lines is shown in Fig. |3l 



trum of Ref. |21|. In addition, the lower Hubbard band is 
broadened due to the opening of the p — d decay chan- 
nel for the d-holes introducing a pronounced asymmetry 
between the upper and lower Hubbard bands ^ (see the 
middle panel of Fig. [2]). Reciprocally, the p bands pos- 
sess a finite k-dependent width due to the coupling to the 
correlated d bands. While full quantitative comparison 
to the ARPES data is not possible because of missing 
dipole matrix elements, the k-dependent broadening of 
the p bands (see the upper panel of Fig. [2]) agrees well 
with the experimental observations In particu- 

lar, we point out the broadening of the otherwise sharp 
p bands observed around the midpoint of the F-X line 
and a considerable smearing of the lowest two p bands 
near the K point (upper right panel of Fig. [2]), which are 
hardly recognizable in the experimental spectra. 

The low energy bands of charge-transfer systems, es- 



Compare to the symmetric upper and lower Hubbard bands of 
eg symmetry in Fig. 3 of Ref. [ISj where p — d hybridization was 
not included. 



FIG. 3: (color online) Detail of the uppermost valence band 
along the F-K (right) and F-X (left) lines. The top panels 
show the 0-p contribution App(k,uj), while Ni-d contribution 
Add(k.,i^) Add(k,i^) of the eg symmetry is shown in the bot- 
tom panels. 



pecially in the case of high temperature cuprate super- 
conductors, have been subject of numerous theoretical 
investigations. This was initiated by Zhang and Rice [llj 
who constructed an effective t — J Hamiltonian for holes 
doped to the copper-oxygen plane and who introduced 
the notion of a bound state between the p-hole and c?-spin 
known as Zhang-Rice singlet. Using a canonical trans- 
formation of the Hubbard model onto the spin-fermion 
model Eroles et al. found a strong k-dependence in 
the orbital composition of the Zhang-Rice band. Bala et 
al. 
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26| applied a generalized spin-fermion model to a 
2D slab of NiO and obtained qualitative agreement with 
the corresponding part of the ARPES spectrum. They 
found very strong k-dependence of the spectral weight 
in the uppermost valence band. To a good approxima- 
tion the latter corresponds to the p-spectral weight of 
the original multi-band Hubbard model 2J| and is thus 
directly comparable to the results of the present study. 
Like Bala et al., we find a vanishing p spectral weight 
in the uppermost valence band along the F-X ((10) in 
Ref. I25I ) line as well as a substantial p contribution in 
the F-K direction ((11) in Ref. ^ shown in Fig. O 
In addition, we find a rather k-independent d contribu- 
tion to the uppermost valence band, which was suggested 
in Ref. 26 and which reflects the local character of the 
Zhang-Rice bound state, a doublet in the case of NiO. 
Unlike model theories, which are restricted to a special 
part of the Hilbert space, the LDA-I-DMFT scheme pro- 
vides a unified picture of all energy scales and avoids 
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any adjustable parameters. The good agreement with 
the experimental data and the low-energy model theory 
found in this work demonstrates the capability of DMFT 
to adequately describe charge-transfer systems and puts 
the earlier model results on a solid foundation. 

In conclusion, by employing dynamical mean-field the- 
ory combined with LDA electronic structure calculation 
we presented the solution to a long-standing problem, the 
computation of the full valence bandstructure of a charge- 
transfer insulator. We obtained a very good agreement 
with the ARPES data of Shen et al. without any 

adjustable parameters. We found well-separated Zhang- 
Rice bands at the top valence manifold with strongly 
k-dependent orbital composition. Our results clearly 
demonstrate the capability of DMFT to treat, upon ex- 
plicit inclusion oip — d hybridization, the late transition- 
metal oxides and charge-transfer systems in general. 
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